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Extended  Abstract 

Field  observations  of  the  ocean’s  forced  stage  response  to  three  hurricanes,  Norbert 
(1984),  Josephine  (1984)  and  Gloria  (1985),  are  analyzed  and  presented  in  a  storm- 
centered  coordinate  system.  The  data  are  simulated  with  two  different  numerical  ocean 
models. 

All  three  hurricanes  had  a  non-dimensional  speed  of  0(1)  and  produced  a  strongly 
rightward  biased  response  of  the  ocean  surface  mixed  layer  (SML)  transport  and  cur¬ 
rent.  This  asymmetry  arises  because  the  wind  stress  vector  rotates  clockwise  on  the 
right  side  of  the  track  and  remains  nearly  parallel  with  the  inertially  rotating  SML 
current  during  most  of  the  hurricane  passage.  The  maximum  layer-averaged  SML  cur¬ 
rents  varied  from  0.8  m  s_1  in  response  to  Josephine,  which  was  a  fairly  weak  hurricane, 
to  1.7  m  s-1  in  response  to  Gloria,  which  was  much  stronger.  In  these  two  cases  the 
current  amplitude  is  set  primarily  by  the  strength  of  the  wind  stress  and  its  efficiency 
of  coupling  with  the  SML  current,  and  the  depth  of  vertical  mixing  of  the  SML  —  all 
of  which  are  local  phenomena  that  can  be  accounted  for  in  a  one-dimensional  model. 
However,  in  the  Norbert  case  where  the  SML  Burger  number  was  appreciable,  ~  1/2, 
the  SML  current  was  also  affected  by  significant  pressure-coupling  with  the  thermo- 
cline  —  an  inherently  non-local  phenomenon.  The  observations  and  the  results  from 
a  three-dimensional  simulation  model  show  appreciable  upwelling  by  inertial  pumping 
and  strong  thermocline-depth  currents,  up  to  0.3  m  s-1,  under  the  trailing  edge  of 
Norbert. 

The  observed  SML  current  has  a  vertical  shear  in  the  direction  of  the  local  wind 
of  about  0.01  s-1  beneath  the  high  stress  regions  of  the  hurricanes.  This  vertical  shear 
causes  the  surface  current  to  be  larger  than  the  layer- averaged  SML  current  described 
above  by  typically  0.2  ms-1,  and  has  been  simulated  with  some  success  by  a  one¬ 
dimensional  turbulence  closure  model. 
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1  Introduction 


The  ocean’s  baroclinic  response  to  a  moving  hurricane  can  be  thought  to  occur 
in  two  largely  distinct  stages  —  a  forced  stage  response  during  the  actual  storm  pas¬ 
sage  and  a  relaxation  stage  response  during  the  following  several  weeks.  These  are 
characterized  by  quite  different  physical  processes. 

The  forced  stage  response  is  a  mainly  local  (depth  and  time  dependent),  non- 
adiabatic  response  of  the  upper  ocean  to  intense  wind  stress.  The  ocean’s  response 
includes  very  large  surface  waves  and  currents,  and  substantial  cooling  of  the  surface 
mixed  layer  (SML)  due  mainly  to  vertical  mixing  (Price,  1981).  The  SML  cooling  can 
be  observed  by  a  variety  of  means,  including  satellite  infrared  imagery,  and  is  now  a 
well- documented  phenomenon  (Black,  1983;  Stramma  et  al.,  1986). 

The  relaxation  stage  response  is  an  inherently  non-local  (three-dimensional  and 
time  dependent),  adiabatic  energy  dispersion  by  near-inertial  frequency  internal  waves 
that  form  a  spreading  wake  behind  a  moving  hurricane  (Geisler,  1970;  Price,  1983). 
Energy  of  the  SML  currents  is  dispersed  vertically  into  the  thermocline,  and  laterally 
away  from  the  storm  track.  These  inertial  waves  are  readily  observed  only  by  direct 
current  observations,  and  thus  only  a  few  case  studies  are  now  available  (Brooks,  1983; 
Shay  and  Elsberry,  1987;  Brink,  1989;  Shay  et  al.,  1989). 

The  connection  between  the  forced  stage  and  the  relaxation  stage  response  can  be 
found  in  the  horizontal  structure  of  the  forced  stage  wind-driven  currents  in  the  SML. 
In  this  paper  we  continue  an  analysis  of  the  forced  stage  response  to  three  open  ocean 
hurricanes  begun  by  Sanford  et  al.  (1987)  (hereafter,  Part  I).  Part  I  described  two  field 
studies  carried  out  in  fall  1984  around  hurricane  Norbert  in  the  eastern  North  Pacific 
Ocean  and  Josephine  in  the  Sargasso  Sea.  In  this  follow-on  report  we  are  able  to  include 
a  third  data  set  acquired  from  a  very  similar  field  study  carried  out  in  fall  1985  with 
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hurricane  Gloria,  also  in  the  Sargasso  Sea.  The  data  collection  and  analysis  methods 
applied  in  the  Gloria  data  were  identical  to  those  described  in  Part  I  (a  brief  review  of 
analysis  methods  and  a  tabulation  of  Gloria  data  are  in  the  Appendix).  These  three 
data  sets  are  unique  in  showing  the  horizontal  and  vertical  structure  of  the  forced  stage 
response,  while  also  providing  the  means  to  estimate  the  two-dimensional  and  time- 
dependent  field  of  surface  wind  stress  that  are  required  for  simulations.  Our  scientific 
goals  are  to  use  these  data  to: 

i)  describe  and  interpret  the  structure  of  the  forced  stage  response,  and 
consider  how  SML  and  thermocline  currents  will  begin  to  evolve  in  the 
relaxation  stage,  and, 

ii)  compare  the  field  data  with  the  results  from  two  numerical  ocean  models 
that  attempt  to  simulate  the  forced  stage,  baroclinic  response. 

The  model  solutions  have  proven  to  be  a  valuable  guide  to  interpreting  the  ob¬ 
servations.  They  help  to  reveal  the  large-scale  pattern  underlying  the  discrete  AXCP 
(air-deployed  expendable  current  profiler)  observations,  and  they  help  also  by  showing 
the  dynamical  relationships  among  variables  —  SML  current  and  sea  surface  cooling, 
and  SML  current  divergence  and  thermocline  current  —  that  would  be  hard  to  appre¬ 
ciate  from  the  data  alone. 

The  practical  goal  of  these  studies  was  to  observe  and  model  the  near  surface  cur¬ 
rents  that  occur  during  the  forced  stage  response.  These  currents  are  an  important 
parameter  in  the  design  of  offshore  structures  that  could  be  subject  to  hurricane  con¬ 
ditions.  As  a  rule,  offshore  structures  are  designed  to  withstand  the  worst  combination 
of  storm-driven  waves  and  currents  that  could  be  expected  during  a  span  of,  typically, 
100  years.  It  is  thus  necessary  to  extrapolate  or  model  the  design  current  from  observed 
cases  which  are  very  likely  to  be  less  severe. 

We  begin  in  Section  2  with  a  description  of  the  hurricane  wind  fields,  and  continue 
in  Section  3  with  the  ocean  initial  condition,  both  of  which  are  important  in  determining 
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the  amplitude  and  structure  of  the  response.  In  Section  4  we  briefly  introduce  two 
numerical  ocean  models  used  to  simulate  the  ocean  response.  The  first  model  is  three- 
dimensional  and  simulates  the  full  horizontal  structure  of  the  response,  but  at  the 
expense  of  representing  vertical  structure  with  only  four  layers  (Price,  1981;  termed  ‘3- 
D’).  The  second  model  has  very  high  resolution  in  the  vertical  dimension  only,  and  is 
useful  for  examining  the  vertical  shear  of  current  within  the  SML  (Mellor  and  Durbin, 
1975;  termed  T-D’).  In  Section  5  we  describe  the  horizontal  structure  of  the  forced 
stage  SML  response,  and  compare  observations  with  the  3-D  model  simulations.  In 
Section  6  we  consider  the  coupling  between  the  SML  and  the  thermocline  that  makes 
up  the  relaxation  stage  response.  The  vertical  structure  of  the  SML  current  and  density 
are  analyzed  in  Section  7  along  with  1-D  model  simulations.  In  the  final  section,  8,  we 
summarize  results  and  suggest  where  future  research  might  be  useful. 


2  The  Hurricanes 

The  NOAA  P3  aircraft  that  carried  out  the  AXCP  surveys  described  in  Part  I 
also  made  measurements  of  the  flight  level  wind  speed,  wind  direction,  and  pressure. 
Flight  level  varied  considerably  between  the  three  cases;  flight  level  was  1500  m  in 
Norbert,  500  m  in  Josephine,  and  considerably  higher  in  Gloria,  3300  m,  because  of 
severe  turbulence.  Because  the  Norbert  and  Josephine  data  were  taken  within  the 
planetary  boundary  layer  (PBL),  we  can  attempt  to  estimate  surface  winds  from  a 
simple  extrapolation  procedure  which  would  not  be  appropriate  in  the  Gloria  case.  In 
the  latter  case  we  use  the  results  from  a  hindcast  simulation  of  the  hurricane  PBL. 
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2.1  Surface  Stress  Estimation  for  Norbert  and  Josephine 

2.1.1  Flight  Level  and  Surface  Winds 

Flight  level  winds  were  sampled  at  very  rapid  intervals  along  each  of  the  radial 
sections  flown  through  the  hurricanes.  These  data  were  subsampled  to  provide  a  wind 
vector  at  roughly  10  km  intervals  (Figure  11  of  Part  I).  These  subsamples  were  taken 
when  the  aircraft  was  in  a  more  or  less  steady  attitude  to  minimize  measurement 
errors.  The  largest  flight  level  winds  observed  in  Norbert  were  about  115  knots  when 
the  central  pressure  was  about  950  mb,  and  the  largest  winds  in  Josephine  were  about 
85  knots  when  the  central  pressure  was  about  974  mb.  Thus  Norbert  was  a  moderately 
intense  hurricane,  while  Josephine  was  somewhat  weaker. 

To  calculate  surface  stress  we  must  extrapolate  the  flight  level  winds  to  some 
standard  height  near  the  surface,  here  taken  to  be  10  m.  This  extrapolation  is  done  by 
a  simple  reduction  of  amplitude  estimated  from  an  empirical  (Bates)  profile  given  by 
Powell  (1980),  and  by  a  rotation  of  the  wind  vector  toward  lower  pressure  estimated 
from  Frank’s  (1977)  data.  Because  flight  level  was  somewhat  higher  in  Norbert,  the 
amplitude  reduction  and  vector  rotation  are  both  slightly  larger  than  were  appropriate 
for  Josephine,  Table  I.  The  uncertainty  in  the  amplitude  reduction  was  estimated  from 
the  error  bars  on  the  original  figures  of  Powell  (1980). 

2.1.2  Fitting  and  Translation  of  Model  Hurricanes 

The  3-D  model  requires  a  two-dimensional  field  of  wind  stress  at  all  time  steps, 
and  hence  these  observed  winds  have  to  be  analyzed  onto  a  regular  grid.  We  chose 
to  fit  these  data  to  a  model  hurricane  (rather  than  interpolate  onto  a  grid)  because 
the  front  left  half  of  Josephine  was  not  sampled,  and,  also  because  a  fitting  procedure 
yields  a  much  more  portable  result.  The  form  of  the  model  hurricane  was  taken  to  be 
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TABLE  I 


Wind  Analysis  and  Hurricane  Model  Fitting 
Extrapolation  to  10  m  height 


Norbert 

Josephine 

Nominal  flight 

Level,  m 

1500 

500 

Amplitude 

Reduction,  percent 

33  ±  7 

27  ±  7 

Rotation,  deg 

16 

10 

Model  Hurricane  Radial  Profiles 


Radius 

Wind  Speed 

Inflow  An 

normalized  by  Rmax 

normalized  by  Umax 

deg 

0. 

0. 

0 

0.4 

0.1 

2 

0.7 

0.5 

4 

0.8 

0.8 

6 

0.95 

0.95 

i 

1.0 

1.0 

7 

1.35 

0.97 

14 

2.7 

0.72 

23 

4.05 

0.54 

24 

5.4 

0.44 

22 

6.75 

0.4 

21 

S.l 

0.36 

21 

10.8 

0.27 

21 

13.5 

0.23 

21 

27.0 

0. 

20 

5 


Norbert 


Best  Fit  Parameters 


Josephine 


Radius  to 


maximum  wind, 

Rma*  i  km 

20  ±  2 

52  ±  3 

Maximum  wind 
speed,  Umax,  m  s-1 

36  ±  2 

29  ±  2 

Percent  variance 
accounted  for 

93 

94 

Rms  error  of 
the  fit,  m  s-1 

6 

6 

6 


the  composite  hurricane  compiled  by  NOAA/NWS  (1979)  for  use  in  design  studies. 
This  model  hurricane  is  specified  by  wind  speed  and  inflow  angle  as  a  discrete  function 
of  radius  in  Table  I.  The  radius  is  normalized  by  the  radius  to  maximum  wind  speed, 
Rmaxi  and  wind  speed  is  normalized  by  the  maximum  wind  speed,  Umax  (dimensional 
values  given  below). 

Translation  of  a  hurricane  will  induce  an  asymmetry  of  wind  speed,  with  larger 
values  occurring  on  the  right  side  of  the  track  where  the  cyclonic  winds  and  the  hur¬ 
ricane  translation  add  constructively.  To  account  for  this,  the  along  track  component 
of  wind  in  the  model  hurricane  is  increased  on  the  right  side  (and  decreased  on  the 
left  side)  by  an  amount  1.17  U fi63  suggested  by  NOAA/NWS  (1979)  (and  ignoring  the 
dimensions  on  Uh  for  this  purpose  only).  Because  Norbert  and  Josephine  were  moving 
fairly  slowly  when  surveyed,  Uh  —  4ms-1,  the  estimated  asymmetry  in  wind  speed 
induced  by  the  translation  is  also  fairly  small,  only  about  2.5  m  s-1,  but  does  make  a 
noticeable  difference  in  estimated  stress. 

The  model  hurricane  (including  the  asymmetry  noted  above)  was  fitted  to  the  10  m 
winds  by  varying  Rmax  and  Umax  to  minimize  the  mean  square  difference  between  the 
observed  winds  and  the  model  hurricane  winds.  The  best  fit  7?max  and  Umax  for  each 
hurricane  are  listed  in  Table  I  along  with  the  percent  variance  accounted  for  by  the  best 
fit  model,  and  the  root  mean  squared  (rms)  wind  velocity  that  could  not  be  accounted 
for  by  the  best  fit.  In  both  cases  the  percent  variance  accounted  for  was  in  excess 
of  90%,  suggesting  that  the  NOAA/NWS  (1979)  model  hurricane  was  appropriate  for 
Norbert  and  Josephine.  Visual  comparison  of  the  observed  and  fitted  winds  gives  the 
same  impression.  There  are,  of  course,  mesoscale  variations  of  the  observed  wind  which 
are  not  included  in  the  model  hurricane.  However,  these  are  small  compared  to  the 
largest  winds,  and  do  not  have  a  systematic,  hurricane-scale  pattern. 
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2.1.3  Calculation  of  Surface  Stress 


The  surface  stress,  r,  was  calculated  from  the  model  10  m  wind,  t/10,  using  the 
usual  bulk  transfer  formula,  t  —  paCdUioUi0,  where  pa  is  the  density  of  air,  and  the 
drag  coefficient  Cd  is  the  Large  and  Pond  (1981)  neutral  form;  Cd  =  1.14  x  10~3  if 
U10  <  10  m  s-1,  and,  Cd  =  (0.49  +  0.065t^o)  x  10~3  if  U\o  >  10  m  s-1.  This  form  is 
very  similar  to  that  inferred  from  hurricane  wind  observations  using  the  ageostrophic 
method  by  Miller  (1964)  (and  see  also  Powell,  1980).  In  these  cases  the  wind  speed 
dependence  of  Cd  is  important  inasmuch  as  the  calculated  Cd  at  the  largest  wind  speeds, 
35  m  s-1,  is  about  2.7  X  10-3,  or  roughly  twice  the  usual  low  wind  speed  value. 

2.1.4  Norbert  and  Josephine  Hurricane  Tracks 

Both  Norbert  and  Josephine  moved  along  rather  complicated  tracks  as  they  passed 
over  the  study  regions.  Model  experiments  have  shown  that  course  or  translation  speed 
changes  are  important  if  they  cause  a  change  in  stress  direction  (as  seen  from  the  ocean) 
of  more  than  about  20  deg,  or  equivalently,  if  they  cause  a  change  in  hurricane  residence 
time  of  more  than  about  one  hour.  In  these  cases  then,  it  was  necessary  to  translate 
the  model  hurricanes  over  the  observed  tracks,  which  were  obtained  from  satellite  and 
aircraft  reconnaissance  by  the  NOAA  P3  (Figure  1).  The  3-D  model  integration  was 
begun  when  the  hurricanes  were  well  outside  the  survey  region,  and  continued  until  the 
hurricanes  reached  the  central  positions  of  the  surveys  given  in  Table  II.  At  that  point 
the  integrations  were  stopped,  and  the  solutions  saved. 


2.2  Surface  Stress  Estimation  for  Gloria 

Because  the  flight  level  in  Gloria  was  at  3300  m  and  above  the  planetary  boundary 
layer,  we  have  not  applied  the  simple  extrapolation  and  model  fitting  procedure 
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TABLE  II 


Position  and  Speed  of  Hurricane  in  Study  Region 


Norbert 

Josephine 

Gloria 

Time  and  date, 

0112,  24  Sep. 

0941,  11  Oct. 

0700,  26  Sep. 

UTC 

Position 

19.41N,  109.08W 

29.41N,  72.13W 

28.75N,  74.98W 

Translation 

4.1  ±0.3 

3.5  ±0.2 

6.8  ±  1.0 

speed,  Uh  ,  m  s-1 

Course,  deg  T 

320  ±5 

10  ±5 

333  ±  5 

10 


described  above.  Instead,  the  Gloria  stress  field  was  estimated  for  this  project  by  Dr. 
V.  J.  Cardone  of  Oceanweather  Inc.  using  a  numerical  model  of  the  planetary  boundary 
layer  below  a  translating,  atmospheric  vortex.  The  time-dependent  pressure  field  of  the 
vortex  and  the  synoptic  scale  environment  were  specified  from  observations  at  30  min 
intervals.  The  model  then  computed  the  winds  and  the  surface  stress  by  means  of  a 
similarity  theory  of  the  marine  boundary  layer,  and  provided  the  field  of  surface  stress 
on  a  regular,  fixed  grid. 

A  series  of  validation  tests  (Cardone  and  Ross,  1979;  Cardone  et  al.,  1980;  For- 
ristall  et  ah,  1977;  Forristall  et  ah,  1978)  indicate  that  this  method  produces  hourly 
average  surface  winds  that  are  generally  accurate  to  within  about  2  m  s_1  in  magni¬ 
tude,  and  20  degrees  in  direction.  This  is  comparable  to  the  uncertainty  inherent  in 
the  extrapolation  procedure  described  above. 

We  have  also  had  the  Norbert  and  Josephine  pressure  data  sets  rim  through  the 
Oceanweather  analysis,  and  found  that  the  stress  fields  were  very  similar  to  those 
estimated  from  the  extrapolation  procedure.  Comparisons  showed  that  the  maximum 
stress  occurred  at  slightly  larger  radius  (by  about  10  to  15%),  but  that  the  maximum 
stress  value  was  almost  identical,  as  was  the  overall  storm  size.  The  end  result  was  that 
the  ocean  models  gave  virtually  the  same  simulations  as  when  driven  with  the  wind 
fields  described  above.  On  that  basis,  we  presume  that  the  Gloria  wind  stress  fields  are 
equivalent  to  those  from  the  analysis  of  Norbert  and  Josephine. 

2.3  Comments  on  the  Wind  Stress 

In  each  model  hurricane  the  largest  stress  occurs  in  the  right  front  quadrant  because 
of  the  asymmetry  induced  by  translation  (Figure  2,  right  panels).  The  left  to  right 
asymmetry  of  stress  is  about  35%,  which  is  enough  to  make  a  detectable  difference  in 
the  simulated  currents.  However,  as  we  will  see  in  Section  5,  the  rightward  asymmetry 
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STRESS.  Po 


DISTANCE  FROM  eye,  km 


Figure  2.  Wind  stress  for  hurricanes  Norbert,  Josephine,  and  Gloria  (top  to  bottom) 
shown  as  friction  velocity  in  plan  view  (left  column)  and  in  cross  section  through  the 
hurricanes  (right  column;  solid  line  is  a  slice  along  the  track,  dashed  line  is  across  the 
track).  In  all  plan  view  figures  the  hurricane  is  centered  in  the  figure,  and  coordinates  are 
rotated  so  that  hurricane  translation  is  in  the  positive  along  track  direction.  The  dotted 
lines  are  the  hurricane  tracks  at  2000  second  intervals. 
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in  the  current  response  is  much  larger  than  this  (roughly  a  factor  of  five)  on  account  of 
an  inherent  asymmetry  in  the  coupling  between  the  surface  stress  and  the  near-inertial 
current  in  the  SML. 

It  is  important  to  keep  in  mind  that  the  wind  stress  estimated  by  either  of  the 
methods  described  above  is  most  directly  related  by  calibration  of  the  aerodynamic 
formula  to  the  stress  in  the  marine  planetary  boundary  layer  (PBL),  and  not  to  the 
stress  within  the  SML  that  we  need  to  know  to  model  currents.  Moreover,  the  bulk 
formula  and  drag  coefficient  are  calibrated  for  more  or  less  steady  conditions  in  which 
surface  waves  can  be  assumed  to  be  a  transparent  intermediary  in  the  momentum 
transfer  from  the  wind  to  the  surface  current  (Stewart,  1974).  In  hurricane  conditions 
surface  waves  are  presumably  quite  non-stationary,  and  could  possibly  lead  to  a  non¬ 
local  momentum  transfer  from  the  PBL  to  ocean  currents.  This  raises  the  following 
question  that  we  can  attempt  to  answer  with  these  data. 

Question  1:  Is  the  hurricane  wind  stress  estimated  from  the  bulk  formula 
absorbed  locally  by  ocean  currents? 

This  will  be  taken  up  in  Section  5  when  we  examine  the  transport  in  the  SML. 

2.4  Comparison  of  the  Hurricanes 

The  hurricanes  can  be  characterized  by  the  three  parameters  Rmax,  Uh ,  and  the 
stress  amplitude,  r,  that  are  listed  in  Table  III  along  with  a  scale  estimate  of  the 
SML  current  amplitude  and  the  upwelling  amplitude  (Price,  1983;  Greatbatch,  1984). 
The  primary  difference  between  these  hurricanes  is  in  their  size;  Norbert  was  a  fairly 
compact  hurricane  having  Rmax  =  20  km,  while  Gloria  was  both  large,  having  Rmax  = 
75  km,  and  also  quite  powerful. 
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TABLE  III 


Hurricane  Parameters 


Norbert 

Josephine 

Gloria 

Uh,  m  s  1 

4.1 

3.5 

6.8 

TLnaxt  km 

20 

52 

75 

r,  Pa  (a) 

4.0 

2.5 

4.0 

/,  s'1 

4.8  x  10“5 

7.1  x  10“5 

7.0  x  10 

hi,  m  (b) 

40 

60 

50 

ft  =  n.  s-  (c) 

0.47 

0.60 

0.86 

(d) 

20 

10 

8 

S  -  — &■— 

2/  Rmax 

2.1 

0.5 

0.6 

M==  +  (e) 

0.5 

0.2 

0.07 

^  ~  Po  h,  Uh  f 

1.5 

0.5 

0.5 

Notes: 

(a)  An  average  of  the  maximum  values  on  the  left  and  right  sides  of  the  hurricane 

(b)  Typical  values  of  z\ 

(c)  This  is  a  scale  for  SML  current.  Even  in  a  linear  case  the  actual  maximum 

current  depends  also  upon  S  and  M,  but  for  S  >  1,  the  maximum  current  is 
about  2U\. 

(d)  Scale  for  the  amplitude  of  inertial  pumping.  Same  comment  holds  as  for  17; 
actual  maximum  is  about  2ijx. 

(e)  The  reduced  gravity  of  the  thermocline,  g',  is  taken  to  be  4  x  10-3<7  in  each 

case. 
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Some  qualitative  aspects  of  the  ocean’s  response  can  be  anticipated  from  the  fol¬ 
lowing  three  non-dimensional  parameters  formed  from  the  characteristic  scales  of  the 
hurricanes  and  the  ocean’s  stratification  (Table  III).  These  are: 

1) 5,  the  non-dimensional  storm  speed,  which  is  the  ratio  of  the  hurricane  residence 
time  to  the  local  inertial  period.  5  is  0(1)  in  each  case,  so  that  the  wind  stress  seen 
from  the  ocean  changes  on  a  time  scale  comparable  to  the  local  inertial  period.  As  a 
consequence,  we  would  expect  that  the  SML  current  should  be  dominated  by  inertial 
motions,  rather  than  geostrophic  currents  which  would  be  more  important  at  smaller  5. 

2)  M,  a  Burger  number  for  the  SML,  which  is  the  ratio  of  the  pressure  gradient 
force  to  the  Coriolis  force  acting  upon  the  SML  current.  Given  the  large  M  in  the  Nor- 
bert  case,  we  would  expect  that  pressure-coupling  between  the  SML  and  thermocline 
would  be  most  pronounced  in  that  case. 

3)  Q ,  a  Rossby  number  for  the  SML  current,  which  is  the  ratio  of  horizontal 
advection  of  momentum  to  the  Coriolis  force.  Again,  the  small  size  of  Norbert  would 
be  expected  to  lead  to  enhanced  non-local  effects  during  the  forced  stage  response. 


3  The  Ocean  Initial  Condition 

In  order  to  analyze  the  data  and  run  the  ocean  models  we  have  to  specify  the  ocean 
initial  condition  over  the  study  regions.  In  principal  this  should  include  the  fields  of 
temperature,  salinity,  and  current.  In  practice  though,  the  aircraft  and  instrument 
resources  required  to  obtain  such  data  were  not  available,  and  even  if  they  had  been,  it 
would  have  been  problematic  to  forecast  the  hurricane  track  well  enough  to  carry  out 
a  useful  survey.  In  the  absence  of  a  complete  initial  survey  we  are  forced  to  make  the 
following  three  assumptions,  the  last  two  of  which  are  dubious  in  the  Josephine  and 
Gloria  cases.  We  assume  that: 
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1)  salinity  played  no  important  role  in  the  density  stratification, 

2)  the  temperature  profile  was  horizontally  homogenous,  and  could  be  estimated 
from  AXCP  temperature  profiles  made  under  the  leading  edge  of  the  hurricanes  (Fig¬ 
ure  3:  Norbert  AXCPs  2,  4,  and  31;  Josephine  AXCPs  20,  21,  7,  3,  and  4;  Gloria 
AXCPs  31,  17  and  18),  and  that 

3)  pre-hurricane  currents  in  the  survey  regions  were  negligible. 

The  estimated  initial  temperature  profiles  are  listed  in  Table  IV,  along  with  un¬ 
certainties  estimated  from  the  variability  among  the  AXCPs  noted.  This  uncertainty 
does  not  include  possible  inhomogeneity.  As  noted  in  Part  I,  satellite  imagery  from  the 
Sargasso  Sea  (NWS/NESS  Oceanographic  Analysis)  shows  significant  pre-hurricane 
horizontal  variability  due  probably  to  the  subtropical  front  (Voorhis,  1969;  and  see 
Black  et  al.,  1988  for  direct  observations  relevant  to  the  Josephine  case).  Evidence  of 
this  in  AXCP  data  is  coolest  surface  temperatures  found  ahead  of  the  hurricanes  (north 
of  the  front),  which  can  not  be  attributed  to  the  hurricane  response,  and  thermocline- 
depth  isotherm  displacements  and  currents  which  are  similarly  incompatible  (coolest 
temperatures  and  largest  thermocline-depth  currents  ahead  of  the  hurricane).  As  a 
result,  the  interpretation  of  the  thermal  field  and  thermocline  current  data  from  the 
Sargasso  Sea  cases  (Josephine  and  Gloria)  is  problematic.  For  those  fields  we  emphasize 
the  Norbert  data  set  which  evidently  had  much  less  pre-hurricane  variability. 

The  neglect  of  initial  currents  will  introduce  an  error  or  uncertainty  in  our  analysis 
about  which  we  can  say  very  little.  We  suspect  that  the  resulting  error  is  as  large  as 
the  uncertainly  associated  with  the  AXCP  measurement  or  analysis  method  (about 
0.2  m  s”1,  Part  I).  One  can  assume  that  SML  currents  of  0(0.2  ms-1)  are  likely  to  be 
found  at  any  open  ocean  site,  and  we  have  to  assume  that  such  currents  could  be 
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Figure  3.  AXCP  profile  numbers. 
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TABLE  IV 


Initial  Temperature  Profiles 


Norbert 

Josephine 

Gloria 

Latitude,  deg. 

19 

29 

29 

Sea  surface 
temperature,  C 

28.5 

27.0 

2S.5 

Mixed-layer 
thickness,  m 

25  ±5 

35  ±  10 

40  ±  5 

Temperature  jump 
at  SML,  C 

1  ±  1 

1  ±  1 

3  ±  1 

Layer  two 
temp,  grad.,  C  m*1 

0.16  ±  0.04 

0.06  ±  0.02 

0.06  ±  0.02 

Layer  two 
thickness,  m 

100 

180 

200 

Layer  three 
temp,  grad.,  C  m_l 

0.016 

0.005 

0.005 

Layer  three 
thickness,  m 

350 

375 

375 

Layer  four 
temp,  grad.,  C  m_l 

0.016 

0.03 

0.02 

Layer  four 
thickness,  m 

380 

450 

580 
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present  and  superimposed  upon  the  total  current  measured  by  the  AXCPs.  The  net 
uncertainty  in  interpreting  the  AXCP  current  data  as  if  it  were  strictly  a  hurricane 
response  is  thus  of  0(0.3  m  s_1).  The  saving  feature  of  this  study  is  that  the  largest 
hurricane-driven  SML  currents  are  0(1  m  s-1),  and  are  large  enough  to  stand  clearly 
above  almost  any  other  background  variability.  Thus  the  pattern  of  the  hurricane- 
driven  SML  current  is  quite  clear  in  all  three  of  our  data  sets.  Moreover,  there  is  a 
significant  difference  in  SML  current  between  the  three  data  sets  that  can  be  readily 
interpreted  as  a  consequence  of  hurricane  size  or  intensity. 


4  Numerical  Models 

Two  quite  different  kinds  of  numerical  ocean  models  are  used  to  simulate  the  re¬ 
sponse  to  the  hurricanes.  The  first  is  a  fully  three-dimensional,  primitive  equation 
model  described  by  Price  (1981)  (3-D),  and  the  second  is  a  one-dimensional  boundary 
layer  model  described  by  Mellor  and  Durbin  (1975)  (1-D).  These  models  complement 
one  another  in  that  the  3-D  model  simulates  the  full  horizontal  structure  of  the  re¬ 
sponse,  but  attempts  to  represent  vertical  structure  with  only  four  layers  of  which  the 
uppermost  is  the  SML.  The  1-D  model  ignores  horizontal  structure  and  thus  non-local 
effects  altogether,  but  attempts  to  simulate  the  detailed  vertical  structure  within  the 
SML  which  the  first  model  ignores. 

As  we  will  show  from  the  data  in  the  following  sections,  the  forced  stage  response 
is  fundamentally  a  local  process,  but  one  which  has  considerable,  interesting  horizon¬ 
tal  structure.  While  non-local  dynamical  effects  are  clearly  present,  especially  in  the 
Norbert  case,  the  non-local  effects  are  generally  not  of  leading  importance  in  the  forced 
stage  response,  so  that  a  1-D  model  can  produce  useful  simulations  of  SML  currents 
as  demonstrated  by  Martin  (1982).  Similarly,  while  there  is  some  vertical  shear  within 
the  SML,  it  is  not  so  large  as  to  complicate  the  interpretation  of  the  3-D  model  results. 
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(And  continuing  this,  we  see  no  evidence  whatever  that  shear  and  non-local  effects 
interact.)  Hence  it  is  quite  possible  and  in  some  ways  efficient  to  consider  separately 
the  questions  of  horizontal  and  vertical  structure  and  dynamics. 

Since  both  of  these  models  have  been  described  in  detail  before,  the  present  dis¬ 
cussion  is  limited  to  a  consideration  of  which  processes  can  be  regarded  as  well  resolved 
and  reliable,  and  which  model  assumptions  and  parameterizations  can  and  should  be 
tested  with  the  available  data. 

4.1  Three-Dimensional  Layered  Model 

Physical  assumptions  and  simplifications  made  in  the  Price  (1981)  model  are  as 
follows: 

1)  Density  is  represented  by  a  linear  state  equation,  and  variation  of  density  is  ig¬ 
nored  except  in  the  hydrostatic  pressure  equation  (Boussinesq  approximation).  Though 
the  model  can  readily  (and  generally  should)  include  salinity  in  the  density  equation, 
in  this  case  no  salinity  data  were  available,  and  salinity  is  ignored.  None  of  these 
assumptions  is  thought  to  be  serious. 

2)  The  Coriolis  parameter  /  is  assumed  constant,  as  is  the  density  stratification 
within  the  thermocline.  These  are  acceptable  for  short  times  after  the  hurricane  passage 
(less  than  a  few  days),  or  before  the  inertial  wave  wake  has  had  time  +o  disperse 
vertically  and  horizontally  away  from  the  track  (Geisler,  1970;  Price,  1983;  Gill,  1984). 

3)  The  sea  surface  is  taken  to  be  a  rigid  lid  to  exclude  fast  barotropic  waves,  and 
the  abyssal  ocean  is  taken  to  be  infinitely  deep  (reduced  gravity  approximation).  These 
approximations  are  appropriate  for  the  deep,  open  ocean  study  sites  considered  here 
where  the  barotropic  current  response  is  expected  to  be  very  small,  0(0.02  ms-1) 
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(Geisler,  1970;  Cooper  and  Thompson,  1989a,  1989b),  compared  to  the  wind-driven 
current  in  the  SML  which  is  0(1  m  s-1). 


4)  As  we  have  noted  already,  vertical  structure  in  this  model  is  represented  by 
only  four  discrete  layers  (Price,  1981,  Figure  4).  The  SML  and  the  transition  layer  are 
represented  by  the  top-most  layer,  which  is  in  direct  contact  with  the  hurricane.  The 
thermocline  is  represented  by  the  lower  three  layers.  Coupling  between  the  SML  and 
the  thermocline  occurs  only  by  means  of  the  hydrostatic,  baroclinic  pressure  gradient 
produced  by  inertial  pumping,  which  is  a  well  resolved  process  in  the  model. 

Question  2:  Can  the  forced  stage  response  and  the  early  relaxation  stage 
response  be  simulated  by  a  model  having  only  a  few  discrete  layers? 

This  is  taken  up  in  Section  6  when  we  examine  the  vertical  structure  of  the  thermocline 
response. 

5)  Finally,  the  only  subgridscale  (unresolved)  process  treated  in  the  model  is  upper 
ocean  vertical  mixing,  which  in  a  layered  model  is  represented  by  SML  entrainment. 
Entrainment  is  a  crucial  process  in  the  forced  stage  response,  causing  significant  changes 
in  SML  thickness,  and  is  the  main  process  causing  sea  surface  cooling  (Price,  1981; 
Martin,  1982;  Cornillon  et  al.,  1987).  In  this  model  the  entrainment  velocity,  We ,  is 
parameterized  by  the  function, 


We  =  5  x  10  4SV  Rv  4  ,  where  (1) 

_  gSphi 
V  PobV’ 

is  the  bulk  Richardson  number,  g  is  the  acceleration  of  gravity,  8p  is  the  density  jump 
across  the  base  of  the  SML,  p0  is  a  reference  density,  and  8V  =  V,  —  V3  is  the  change 
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of  current  across  the  base  of  the  SML  (layer  2  is  the  upper  thermocline  layer  in  the 
3-D  model).  This  form  gives  significant  entrainment  rates  only  when  Rv  is  less  than 
1,  which  in  practice  occurs  when  the  SML  current,  V,,  is  strongly  accelerated  by  the 
wind  stress. 

Question  3:  Is  the  SML  bulk  Richardson  number  less  than  1  in  regions  of 
strong  entrainment? 

We  will  take  this  up  in  Section  5.3  when  we  examine  SML  cooling  and  the  Richardson 
number  estimated  from  the  AXCP  data. 

The  3-D  model  was  implemented  on  a  grid  with  horizontal  resolution  of  15  km,  and 
integrated  with  a  time  step  of  900  sec,  both  of  which  are  sufficient  to  yield  a  converged 
solution.  There  is  no  explicit  diffusivity  in  the  model,  and  numerical  diffusivity  is 
negligible  for  the  short  duration  of  these  simulations. 


4.2  One-Dimensional  Model 

The  1-D  boundary  layer  model  used  here  was  developed  by  Mellor  and  Durbin 
(1975)  to  simulate  the  wind-forced  evolution  of  upper  ocean  current  and  temperature 
profiles. 

Physical  assumptions  and  simplifications  sire  as  follows. 

1)  Turbulent  fluxes  of  momentum  and  heat  within  the  SML  and  transition  layer 
are  parameterized  by  an  eddy  viscosity  and  an  eddy  diffusivity.  Hence,  this  model 
can  simulate  vertical  shear  within  the  SML.  In  this  particular  version  of  the  model, 
these  exchange  coefficients  are  calculated  from  relatively  simple  algebraic  equations, 
making  this  a  Level  2  scheme  in  the  terminology  of  Mellor  and  Yamada  (1982).  The 
simplifications  in  a  Level  2  scheme  include  the  neglect  of  turbulent  vertical  advection 
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and  diffusion  of  turbulent  energy.  These  terms  should  be  relatively  small  in  a  wind 
driven  boundary  layer  in  which  the  turbulence  production  is  dominated  by  shear. 

2)  The  model  is  one  dimensional,  and  thus  ignores  all  of  the  non-local  dynamics 
that  make  up  the  relaxation  stage  response  —  inertial  pumping,  horizontal  pressure 
gradients,  and  horizontal  advection.  In  the  extreme  cases  where  these  effects  can  be 
fairly  large  (small  intense  hurricanes,  such  as  Norbert),  this  model  alone  might  not 
yield  accurate  simulations  of  the  SML  current.  However,  these  non-local  effects  gen¬ 
erally  have  a  large  vertical  scale  and  contribute  a  nearly  depth-independent  current 
perturbation  to  the  directly  wind-driven  SML  current.  Hence,  it  is  plausible  to  treat 
separately  the  non-local  effects  and  the  issue  of  vertical  shear  within  the  SML. 

The  eddy  coefficients  are  the  product  of  the  turbulent  kinetic  energy,  a  length 
scale  calculated  from  the  first  vertical  moment  of  the  turbulent  energy,  and  a  stability 
function  dependent  upon  the  gradient  Richardson  number.  The  strong  dependence 
upon  a  gradient  Richardson  number  causes  this  model  to  have  a  mixing  response  to 
wind-forcing  that  is  in  some  respects  very  similar  to  that  of  the  3-D  model  in  which 
the  mixing  depends  upon  a  bulk  Richardson  number.  That  is,  this  model  also  gives 
strong  vertical  mixing  in  circumstances  where  the  wind  stress  accelerates  the  upper 
layer  current  to  large  values  and  thereby  produces  low  gradient  Richardson  numbers 
over  a  thick,  upper  ocean  layer. 

Given  the  eddy  viscosity  and  diffusivity  as  a  function  of  depth,  the  one-dimensional 
equations  for  turbulent  transport  can  then  be  solved  numerically.  In  turn,  the  eddy 
coefficients  can  be  calculated  from  the  algebraic  equations.  The  system  is  formally 
closed,  but  iteration  is  required  in  order  to  find  eddy  coefficients  which  are  consistent 
with  the  mean  flow.  The  high  wind  stress  imposed  by  the  hurricanes  caused  stability 
problems  in  the  iterative  solution  until  two  smoothing  operations  were  introduced. 
First  the  updated  eddy  parameters  were  smoothed  with  a  three  point  operator,  and 
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then  the  average  of  the  new  and  old  values  of  the  eddy  parameters  was  used  in  the  bulk 
flow  computations.  Usually,  fewer  than  ten  iterations  were  required  for  convergence. 
In  these  simulations  we  used  a  vertical  grid  size  of  4  m  and  a  time  step  of  1800  sec. 

The  Mellor  and  Durbin  (1975)  model  has  been  applied  extensively,  and  verified 
against  several  upper  ocean  data  sets,  including  current  meter  measurements  in  Hur¬ 
ricane  Eloise  (Martin,  1982).  Our  AXCP  observations  are  particularly  appropriate  for 
testing  this  model  since  they  give  detailed  vertical  profiles  of  current  and  temperature. 
A  specific  question  that  we  address  here  is 

Question  Does  the  1-D  model  give  realistic  solutions  for  vertical  shear 
of  current  within  the  SML  and  transition  layer  ? 

This  will  be  taken  up  in  Section  7. 


5  Horizontal  Structure  of  the  Surface  Mixed-Layer 
Response 

In  this  section  we  describe  the  horizontal  structure  of  the  response  within  the  SML, 
and  compare  the  3-D  model  solutions  to  the  observations.  We  begin  with  a  discussion 
of  the  wind-driven  transport  in  the  SML  because  it  has  simpler  dynamics  than  does 
the  SML  current,  and  because  the  observed  transport  can  be  used  to  make  a  fairly 
sensitive  test  of  the  estimated  wind  stress. 


5.1  Wind-Driven  Transport 

The  wind-driven  (volume)  transport  evaluated  from  the  three-layer  model  fit  to 
AXCP  observations  is  estimated  to  be  D  =  V}hi  +  V2/i2  since  layer  2  is  the  transition 
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layer  over  which  the  SML  current  matches  to  the  pressure-driven  current  in  the  thermo- 
cline.  The  current  in  layer  2  is  nearly  always  found  to  be  parallel  to  the  SML  current, 
and  thus  the  wind-driven  transport  is  very  nearly  parallel  to  the  SML  current.  The 
transport  evaluated  from  the  3-D  model  solution  is  just  D  =  V ih\,  where  subscript  1 
denotes  the  SML. 

The  transport  vectors  estimated  from  ACXP  data  are  plotted  in  Figure  4  (left)  in  a 
storm-centered,  quasi-synoptic  coordinate  system  in  which  the  hurricane  is  translating 
directly  up  the  page  (details  are  given  in  the  last  section  of  the  Appendix).  Simulations 
from  the  3-D  model  are  plotted  in  exactly  the  same  way  in  Figure  4  (right).  To  the 
extent  that  the  translation  is  steady,  then  the  along-track  coordinate  is  time-like.  The 
major  advantage  of  this  storm- centered  coordinate  system  is  that  it  helps  us  to  compare 
the  three  cases. 

5.1.1  Rightward  Bias 

The  dominant  pattern  in  each  case  is  a  clockwise  turning  and  acceleration  of 
the  transport  vectors  underneath  the  hurricanes,  and  a  striking  rightward  bias  in  the 
amplitude  (noted  also  in  Part  I  and  by  Church  et  al.,  1989).  For  example,  at  a  position 
75  km  to  the  right  of  the  center  of  Norbert  the  observed  transport  is  about  60  m2  s_1, 
while  at  the  same  distance  to  the  left  of  the  track  the  amplitude  is  only  about  5  m2  s-1. 
A  comparable  rightward  bias  occurs  in  the  model  solutions  because  of  an  inherent 
asymmetry  in  the  coupling  between  the  wind  stress  of  a  moving  hurricane  and  the 
wind-driven  SML  transport.  The  transport  (and  the  SML  current)  tend  to  rotate 
inertially  (clockwise  in  this  northern  hemisphere  case).  On  the  right  side  of  the  track 
the  wind  stress  also  turns  clockwise  with  time  when  viewed  from  the  ocean.  When 
5  is  0(1)  as  it  is  here  and  for  mos^  hurricanes,  the  wind  stress  rotation  rate  roughly 
matches  the  rotation  rate  of  the  transport  (Chang  and  Anthes,  1978;  Price,  1981),  and 
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Figure  4.  Transport  within  the  SML  as  estimated  from  AXCP  data  (left  panels),  and  as 
simulated  by  the  3-D  model  (right  panels).  A  transparent  overlay  greatly  facilitates  the 
comparison  of  the  observed  and  simulated  fields. 
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the  end  result  is  that  the  transport  and  wind  stress  remain  roughly  aligned  throughout 
most  of  the  hurricane  passage.  Thus  the  transport  is  accelerated  to  almost  the  largest 
possible  amplitude  on  the  right  side  of  the  track. 

On  the  left  side  of  the  track  the  coupling  between  the  SML  transport  and  the  wind 
stress  is  much  less  efficient  because  the  wind  stress  rotates  anti-clockwise  during  the 
hurricane  passage.  Under  the  leading  edge  of  the  hurricane  the  wind  stress  begins  to 
accelerate  the  resting  ocean,  but  within  about  four  hours  the  inertial  rotation  of  the 
transport  and  the  anti-clockwise  rotation  of  wind  stress  cause  the  current  and  wind 
stress  to  be  roughly  anti-parallel  during  most  of  the  hurricane  passage.  The  result  is 
that  transport  never  becomes  very  large  on  the  left  side  of  the  track. 

This  left-to-right  asymmetry  in  transport  amplitude  is  enhanced  by  the  stronger 
wind  stress  amplitude  that  also  occurs  on  the  right  side  of  a  moving  hurricane  (Sec¬ 
tion  2.3).  However,  the  factor  of  five  difference  in  transport  amplitude  across  the 
hurricane  track  is  due  overwhelmingly  to  the  asymmetric  rotation  of  the  wind  stress 
noted  above  rather  than  the  comparatively  small  asymmetry  of  wind  stress  magnitude. 

Notice  that  the  very  simple  dynamics  that  lead  to  the  left-to-right  asymmetry  in 
transport  involve  only  the  local  wind  stress  and  the  local  current  or  transport.  Thus 
the  1-D  model  can  also  simulate  this  asymmetry,  if  driven  with  the  appropriate  time- 
varying  wind  stress  (noted  also  by  Martin,  1982),  and  indeed  the  1-D  model  gives  a 
roughly  comparable  result  for  simulated  SML  currents  and  transports  (more  on  this  in 
Section  7). 

The  storm- centered  coordinate  system  also  helps  to  reveal  the  case-to-case  dif¬ 
ferences  in  the  horizontal  patterns  of  transport,  most  of  which  can  be  attributed  to 
variations  in  the  hurricane  size  or  track.  Notice  particularly  that  the  overall  width 
of  the  region  of  strong  response  is  considerably  bigger  in  the  Gloria  case  than  in  the 
Norbert  case,  and  not  surprisingly  is  directly  proportional  to  the  size  of  the  hurricanes. 
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The  Norbert  case  also  shows  a  region  of  strong  response  to  the  left  rear  which  is  not 
seen  in  the  other  cases.  This  is  due  to  the  cyclonically  curving  track  that  Norbert  took 
upon  entering  the  survey  region  (Figure  1). 

5.1.2  A  Test  of  the  Estimated  Wind  Stress 

For  a  hurricane  of  a  given  size  and  moving  along  a  given  track,  the  simulated 
transport  is  almost  directly  proportional  to  the  wind  stress  amplitude;  it  is  almost 
completely  independent  of  model-specific  physics  and  parameterizations,  and  in  par¬ 
ticular  the  transport  amplitude  does  not  depend  upon  vertical  mixing  (entrainment). 
Simulated  transport  is  affected  somewhat  by  the  pressure-coupling  between  the  SML 
and  the  thermocline,  but  as  noted  in  the  discussion  of  models  in  Section  3,  this  pressure 
coupling  is  a  fairly  well  resolved  process  in  the  3-D  model.  Hence,  we  would  expect 
any  other  three-dimensional,  primitive  equation  model  driven  by  the  same  hurricanes 
to  yield  virtually  the  same  solutions  for  wind-driven  transport  that  we  show  here  from 
our  3-D  model.  A  comparison  of  the  observed  and  simulated  transport  thus  makes  a 
fairly  sensitive  test  of  the  wind  stress  in  near  isolation  from  model-specific  assumptions 
regarding  physics  or  parameterizations. 

To  quantify  the  comparison  of  the  simulated  and  observed  data  we  have  computed 
some  simple  statistics  on  the  observed  and  simulated  transports,  and  on  the  differences 
between  the  two,  Table  V.  For  example,  a  measure  of  the  average  vector  difference 
between  the  observed  and  the  simulated  transport  is  just 


rms  transport  difference 


\ 


N 


i  Y.(Do>-0Di)\ 

;V  1=1 


where  Doi  is  the  observed  transport  vector  at  AXCP  i,  Di  is  the  simulated  transport 
vector  at  the  position  of  AXCP  i  found  by  linearly  interpolating  the  3-D  model  solution, 
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TABLE  V 


Statistical  Comparison  of  Observed  and  Simulated  Velocity  and  Transport 

Norbert 

Josephine 

Gloria 

rms  observed 

37 

24 

68 

transport,  m2  s-1 

(21)* 

rms  transport 

17 

19 

37 

difference,  m2  s'1 

(10) 

average  transport 

7 

0 

-7 

mag.  difference,  m2  s-1 

(-2) 

percent  variance 

80 

38 

70 

of  transport 

(78) 

rms  observed 

0.78 

0.35 

1.05 

velocity,  m  s_1 

(0.32) 

rms  velocity 

0.28 

0.30 

0.44 

difference,  m  s-1 

(0.16) 

average  speed 

0.16 

-0.07 

0.06 

difference,  ms-1 

(-0.07) 

percent  variance 

87 

27 

82 

of  velocity 

(75) 

*  Values  in  parentheses  axe  omitting  AXCPs  2,  3  and  4 
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N  is  the  number  of  AXCPs  in  the  ensemble  average,  and  fi  —  1  here.  From  a  similar 
statistic  computed  over  the  observed  transport  separately  we  can  compute  the  percent 
variance  accounted  for  by  the  simulated  transport  as  100  x  (1  —  ( ^ 'ob^fved 
A  measure  of  the  transport  difference  due  solely  to  differences  in  magnitude  is  just 

1  N 

average  transport  mag.  difference  =  —  ^(|  Doi  |  —  |  D,  I)- 

**  t=i 

These  statistics  show  that  in  the  Norbert  and  Gloria  cases  the  simulated  transport 
can  account  for  about  75%  of  the  variance  in  the  observed  transport  (this  is  true  in  the 
Josephine  case  only  if  we  are  allowed  some  selective  averaging,  more  on  this  below), 
and  that  on  average  over  all  three  cases  the  simulated  transport  has  about  the  same 
magnitude  as  the  observed  transport.  That  is,  there  is  no  evidence  of  a  substantial  bias 
error  in  stress  amplitude,  and  specifically  there  is  no  simple  scaling  of  the  simulated 
transport  fields  by  the  factor  that  succeeds  in  producing  a  smaller  rms  difference  than 
is  obtained  straightaway  with  /?  =  1.  When  the  simulated  transports  are  multiplied  by 
a  scale  factor  0  >  1.2,  or  0  <  0.8,  there  is  a  marked  decrease  in  the  percent  variance 
accounted  for,  and  on  this  basis  we  conclude  that  the  hurricane  stress  fields  used  in 
these  simulations  are  consistent  with  the  observed  ocean  transport  to  within  about 
20%. 

It  remains  that  there  could  be  a  distortion  of  the  stress  field  that  is  largely  cancelled 
by  averaging  over  the  ensemble  of  observations.  Indeed  there  is  an  intriguing  hint  that 
the  transport  difference  may  not  be  completely  random  in  that  the  only  large  errors 
seem  to  occur  in  the  left  rear  quadrant  of  the  hurricane.  By  overlaying  the  observed 
and  simulated  transports  one  can  see  that  the  largest  difference  is  found  at  the  position 
of  Norbert  AXCPs  3,  20  and  21,  and  Josephine  AXCPs  2,  3,  and  4  (Figure  3).  In  the 
latter  case,  satellite  imagery  shows  that  pre-hurricane  thermal  variability  was  present 
and  may  have  had  associated  currents  that  confound  our  attempt  to  interpret  the  total 
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signal  as  if  it  were  hurricane-driven  (the  rationale  for  omitting  these  AXCPs  from  the 
statistics  in  Table  V).  There  is  no  similar  known  feature  in  the  Norbert  case,  however, 
and  so  this  might  also  be  the  result  of  a  sub-hurricane  scale  stress  perturbation  not 
accounted  for  in  these  simulations.  Whether  this  is  due  to  non-stationary  surface  wave 
effects  or  to  some  unresolved  variation  in  the  hurricanes  wind  fields  is  something  that 
we  can  not  tell  from  these  data  alone. 

5.2  SML  Current 

The  field  of  the  layer-averaged  SML  current  (Figure  5)  looks  very  much  like  the 
transport,  and  indeed  all  of  the  previous  discussion  of  horizontal  structure  applies 
equally  well  here.  As  we  implied  in  that  discussion,  the  SML  current  is  a  near-inertial 
motion  having  a  small  blue  shift  of  frequency.  The  clockwise  turning  of  the  current  is 
most  apparent  in  the  Gloria  case  where  the  AXCP  data  span  the  greatest  distance  along 
the  track.  The  wavelength  along  the  track  is  the  inertial  wavelength,  Uh  x  IP(  1  —  v), 
where  IP  is  just  the  local  inertial  period,  and  v  is  the  blue  shift  of  frequency,  typically 
~  0.1  (Price,  1983).  These  data  sets  are  some  of  the  few  ever  made  that  show  the 
horizontal  structure  of  a  near-inertial  current  (see  also  D’Asaro  (1989)  and  Pollard 
(1980)). 

The  SML  current  simulation  by  the  3-D  model  appears  fairly  realistic  in  each  of 
the  cases,  and  certainly  the  large-scale  structure  of  the  current  field  is  reproduced  quite 
well.  The  reasonably  good  comparison  between  the  simulated  and  observed  currents 
encourages  us  to  think  that  the  simulated  current  field  can  be  used  as  an  interpolator 
between  the  discrete  AXCP  data  points.  It  would  appear  that  the  AXCP  sampling, 
though  less  than  we  had  intended,  nevertheless  does  succeed  in  defining  the  overall 
structure  of  the  forced  stage  response  of  the  SML  current. 
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Figure  5.  Current  vectors  within  the  SML.  Note  that  the  qualitative  pattern  is  very 
similar  to  that  of  he  transport,  and  that  the  largest  currents  are  to  the  right  of  the  track 
by  about  25  km  (Norbert)  to  100  km  (Gloria). 
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A  close  look  at  the  current  field  shows  that  it  differs  from  the  transport  in  having 
not  quite  such  a  strong  rightward  bias  in  amplitude.  This  reduced  asymmetry  arises  in 
the  3-D  model  solution  because  enhanced  mixing  on  the  right  side  of  the  track  causes 
the  SML  depth  to  be  somewhat  greater  on  the  right  side  of  the  track,  and  because,  for 
a  given  transport,  the  current  amplitude  is  inversely  proportional  to  SML  depth.  As 
we  will  see  in  the  next  subsection,  there  is  also  evidence  directly  from  the  data  that 
upper  ocean  mixing  is  indeed  enhanced  on  the  right  side  of  the  track. 

The  statistical  comparison  between  the  observed  and  the  simulated  current  (Ta¬ 
ble  IV)  is  about  as  before  for  the  transport.  The  3-D  model  can  account  for  about 
80%  of  the  variance  in  the  observed  layer- averaged  SML  currents.  The  rms  difference 
between  the  simulated  and  observed  currents  is  approximately  0.35  m  s_l,  or  not  much 
larger  than  expected  from  the  total  uncertainty  on  the  observations  (about  0.3  m  s_1 
overall)  (see  remarks  in  Section  3). 

Because  one  goal  of  this  work  was  to  observe  the  maximum  currents  produced  by 
the  hurricanes,  the  flight  plans  included  somewhat  heavier  AXCP  sampling  over  the 
right  side  of  the  tracks.  A  qualitative  comparison  between  the  observed  and  simulated 
fields  suggests  that  the  AXCP  sampling  probably  did  suffice  to  observe  the  strongest 
or  very  nearly  the  strongest  currents  beneath  each  of  these  hurricanes.  The  maximum 
observed  current  varied  by  a  factor  of  almost  2  over  these  three  cases.  In  the  Norbert 
case  the  observed  (simulated)  maximum  current  was  1.10(l.ll)ms_1,in  Josephine  it 
was  0.77  (0.79)  m  s-1,  and  in  Gloria  it  was  1.70  (1.38)  m  s-1.  The  maximum  current 
increases  with  both  the  intensity  and  size  of  the  hurricane  (and  model  studies  suggest 
dependence  upon  other  external  parameters  as  well  (Price,  1981;  Greatbatch,  1984)). 
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5.3  Upper  Ocean  Mixing 


1 


While  we  argued  in  Section  4  that  the  transport  is  nearly  model-independent,  the 
same  is  not  true  for  the  SML  current,  the  SML  depth  and  SML  cooling.  Each  of  these 
is  affected  significantly  by  turbulent  vertical  mixing,  which  has  to  be  parameterized 
in  any  model.  If  vertical  structure  is  represented  by  layers,  as  in  our  3-D  model,  then 
vertical  mixing  may  be  represented  by  SML  entrainment  ,  which  has  been  parameterized 
in  a  variety  of  ways  (e.g.,  the  hurricane  response  models  of  Elsberry  et  al.  (1976)  and 
Chang  and  Anthes  (1978)  give  plausible  results  using  a  form  which  is  quite  different 
from  the  one  used  here).  The  AXCP  data  offer  a  chance  to  make  at  least  a  consistency 
check  on  the  mixing  parameterization  of  the  3-D  model. 

5.3.1  SML  Depth 

When  comparing  the  layered  model  with  oceanic  observations,  we  have  to  keep  in 
mind  that  there  is  a  thick  transition  layer  which  is  a  part  of  the  directly  wind-driven 
layer  extending  well  below  the  literal  mixed  layer  into  the  stratified  fluid  below.  In  the 
AXCP  analysis  the  transition  layer  is  represented  by  layer  2  (Appendix).  Accordingly, 
when  evaluated  from  the  AXCP  data,  an  equivalent  (to  a  layered  model)  SML  depth 
(or  thickness  of  the  wind-driven  layer)  is  estimated  to  be  h\  4-  h2/2,  where  h2  is  the 
transition  layer  thickness  (Figure  6). 

There  is  an  appreciable  small  scale  variability  in  the  estimates  of  the  SML  depth; 
for  example  a  16  m  difference  between  neighboring  AXCPs  3  and  21  from  Norbert. 
This  could  be  due  to  real  horizontal  variations  associated  with  internal  waves,  or  it 
could  also  be  a  result  of  uncertainty  in  estimating  the  transition  layer  thickness  in  the 
AXCP  profiles.  In  either  event,  the  estimated,  observed  SML  depth  is  a  somewhat 
noisy  variable.  There  is  a  fairly  clear  cut  trend  for  SML  depth  to  increase  from  fron  to 
rear.  The  deepest  observed  SML  value  was  62  m  at  AXCP  13  to  the  right  rear  of  the 


34 


ALONG-TRACK  DISTANCE,  km  ALONG- TRACK  DISTANCE,  km  ALONG- TRACK  DISTANCE,  km 

-240  -160  -80  0  80  160  240  -240  -160  -80  0  80  160  240  -240  160  80  0  80  16 


ACROSS-TRACK  DISTANCE,  km 


ACROSS-TRACK  DISTANCE,  km 


Figure  6.  SML  depth  (top  two  panels),  cooling  (middle  two  panels)  and  bulk  Richardson 
number  (bottom  most  panels)  for  the  Norbert  case.  Note  that  the  largest  cooling  found  in 
the  simulation  occurs  to  the  rear  of  the  hurricane,  and  between  two  AXCP  sections. 
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confer.  At  the  same  point  the  3-D  simulation  gives  a  depth  of  about  70  m,  or  somewhat 
deeper.  Within  the  region  of  greatest  mixing  to  the  right  of  the  track,  the  simulated 
SML  depth  roughly  doubled  during  the  passage  of  Norbert. 

5.3.2  SML  Cooling 

The  SML  cooling  (or  SST  cooling)  caused  by  a  hurricane  is  due  almost  entirely  to 
vertical  mixing  (Price,  1981;  Cornillon  et  al.,  1987),  and  should  be  a  particularly  sen¬ 
sitive  diagnostic  of  vertical  mixing.  To  evaluate  the  cooling  from  this  kind  of  one-time 
survey  data,  we  have  no  recourse  but  to  subtract  away  an  estimate  of  the  presumably 
uniform  pre-hurricane  surface  temperature  evaluated  as  part  of  the  initial  condition 
(Section  3,  Table  IV).  The  so-called  cooling  evaluated  in  this  way  is  perhaps  better  de¬ 
scribed  as  a  temperature  anomaly,  and  we  noted  already  in  Section  3  that  it  gives  a  very 
muddled  picture  of  the  cooling  response  to  Josephine  and  Gloria.1  Thus  we  limit  this 
discussion  to  the  Norbert  data  set,  which  apparently  had  much  weaker  pre-hurricane 
variability. 

In  the  Norbert  case  there  is  at  least  a  clear  trend  for  cooling  to  increase  from  front 
to  rear  (Figure  6),  which  must  obtain  if  the  cooling  is  due  to  the  hurricane.  It  appears 
that  there  is  some  rightward  asymmetry  of  cooling  with  the  greatest  values  of  about 
2.2  C  to  the  right  of  the  track.  The  simulation  gives  roughly  comparable  values  at  the 

Cornillon  et  al.  (1987)  have  reported  satellite  infrared  images  of  the  western  North  Atlantic  that 
provide  a  remarkable  view  of  the  sea  surface  cooling  caused  by  hurricane  Gloria.  By  differencing  pre- 
and  post-hurricane  images,  they  were  able  to  isolate  clearly  the  SST  cooling  effect  of  the  hurricane. 
Gloria  caused  pronounced  cooling  along  a  track  extending  from  the  subtropics  all  the  way  north  to  New 
England.  Over  the  Sargasso  Sea,  the  maximum  cooling  was  found  about  100  km  to  the  right  of  the  track, 
or  roughly  where  we  find  the  maximum  SML  current  and  maximum  cooling  response  in  our  models  (1-D 
or  3-D)  (Figures  4  and  5). 
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AXCP  positions,  but  also  suggests  that  larger  cooling  might  have  taken  place  behind 
the  hurricane  where  observations  were  not  taken. 

Overall,  these  comparisons  of  SML  depth  and  cooling  are  rather  unconvincing, 
partly  because  we  can  see  a  reasonably  clear  mixing  signal  in  only  one  of  three  cases, 
but  partly  too  because  the  AXCP  sampling  for  these  variables  appears  not  to  have 
been  as  apt  as  it  was  for  SML  current  and  transport. 

5.3.3  Bulk  Richardson  Number 

In  both  models  the  mixing  effect  of  wind-driven  currents  is  parameterized  by  a 
function  of  a  bulk  or  a  gradient  Richardson  number  (3-D  or  1-D  respectively)  that 
predicts  steeply  increasing  mixing  as  the  Richardson  number  is  forced  below  a  quasi- 
‘criticaT  value.  The  idea  of  a  critical  gradient  Richardson  number  control  on  mixing  is 
well  established  (e.g.,  Mellor  and  Yamada,  1982);  the  equivalent  for  the  bulk  Richardson 
number  is  much  less  so.  The  parameterization  used  here,  Eqn.  (1),  gives  significant 
entrainment  when  Rv  is  less  than  1,  and  very  small  entrainment  for  Ry  greater  than 
1.  If  this  were  indeed  appropriate,  then  we  should  observe  that  Rv  is  somewhat  less 
than  1  in  regions  of  strong  mixing,  and  that  Rv  should  always  be  larger  than  1  where 
mixing  is  very  weak  or  vanishing. 

To  evaluate  Rv  from  the  AXCP  data  we  estimate  the  density  difference  to  be 
8p  =  aST,  where  8T  is  the  temperature  change  across  the  transition  layer  (listed  in 
Table  1  of  the  Appendix),  a  is  the  thermal  expansion  coefficient  (—0.32  kg  m-3  C-1),  the 
SML  thickness  is  estimated  as  noted  above,  and  the  current  difference  is  estimated  to  be 
8V  =  V1  —  Vj,  where  V3  is  the  thermocline  current.  In  profiles  where  the  SML  current 
is  large  (greater  than  about  0.5  m  s-1)  there  is  generally  a  well-defined  transition  layer 
with  a  clear  8T.  In  such  profiles  Rv  is  then  well  defined  from  the  field  data.  However,  at 
locations  where  the  SML  current  is  small  (left  of  the  track  and  ahead  of  the  eye)  there 
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is  only  a  thin  and  sometimes  indistinct  transition  layer,  and  consequent  uncertainty 
in  estimating  ST.  As  a  result,  the  Ry  evaluated  from  the  AXCP  data  is  then  rather 
uncertain,  and  we  suspect  probably  overestimated  in  regions  where  it  would  have  been 
large  (compared  to  1)  anyway. 

The  smallest  Rv  are  found  to  the  right  of  the  track  and  just  behind  the  eye  where 
entrainment  is  expected  to  be  strongest  (Figure  6).  Typical  values  in  this  region  are 
in  the  range  0.5  to  1.0,  which  are  not  inconsistent  with  the  expectations  from  the 
entrainment  parameterization.  Larger  values,  which  are  in  some  cases  very  uncertain 
owing  to  uncertainty  in  the  temperature  change  and  depth  limits  of  the  transition 
layer,  are  found  in  regions  to  the  left  of  the  track,  and  ahead  of  the  hurricane  where 
entrainment  is  weak  or  nearly  vanishing.  Indeed,  all  of  the  Ry  estimates  in  these 
outlying  regions  are  large  compared  to  1.  Thus  the  observed,  bulk  Richardson  number 
has  a  pattern  underneath  the  hurricane  which  is  at  least  roughly  consistent  with  the 
entrainment  parameterization  used  in  the  3-D  model. 

A  similar  pattern  obtains  also  for  the  Josephine  and  Gloria  cases.  An  intriguing 
difference  in  detail  is  that  the  smallest  values  of  Ry  are  slightly  larger  in  the  Josephine 
case,  about  1.0,  and  are  slightly  smaller  in  the  Gloria  case,  as  low  as  0.5.  The  3-D 
model  on  the  other  hand,  tends  to  give  somewhat  more  evenly  distributed  small  values. 
This  may  be  evidence  that  while  a  form  something  like  our  Eqn.  (1)  is  appropriate  for 
these  kinds  of  problems,  there  may  not  be  quite  such  a  steep  roll-off  with  increasing 
Ry  as  suggested  by  (1)  (i.e.,  the  exponent  may  be  less  than  4). 


6  Thermocline  Currents 

An  important  and  somewhat  surprising  result  of  this  study  is  the  observation  of 
substantial  thermodine-depth  currents  under  the  rear  half  of  hurricane  Norbert,  Fig- 
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ure  7.  This  is  evidence  of  an  unusually  strong  pressure-coupling  between  the  SML  and 
thermocline  in  that  case,  and  consistent  with  the  comparatively  large  SML  Burger  num¬ 
ber,  M  =  0.5.  These  large  thermocline  currents  show  that  in  the  Norbert  case  the  forced 
stage  response  and  the  relaxation  stage  response  overlapped  in  time,  and  that  non-local 
dynamics  were  likely  to  be  important  during  the  forced  stage  response.  (There  is  a 
somewhat  weaker  coupling  in  the  other  two  cases  as  well,  but  the  thermocline-induced 
currents  are  apparently  obscured  by  ambient  currents;  evidence  of  this  is  stronger 
thermocline-depth  currents  under  the  leading  edge  of  hurricane  Gloria  than  behind  it.) 

The  3-D  model  gives  a  plausible  simulation  of  the  thermocline  current,  and  helps 
to  show  the  mechanism  of  pressure-coupling  via  inertial  pumping  (see  also  Shay  et  al., 
1989,  who  model  the  deep  response  by  a  summation  of  normal  modes).  With  some 
guidance  from  the  simulation  it  is  easy  to  see  that  the  SML  current  field  (Figure  5)  is 
strongly  divergent  in  a  region  just  behind  the  hurricane;  the  current  to  the  right  of  the 
track  flows  to  the  right  with  a  large  amplitude,  while  the  current  to  the  left  of  the  track 
flows  away  from  the  track  but  is  much  weaker.  This  divergence  acts  to  compress  the 
SML  and  to  upwell  the  thermocline.  Because  the  SML  current  oscillates  with  a  near 
inertial  period,  so  too  does  the  divergence  and  the  associated  up-  and  downwelling, 
which  has  been  termed  ‘inertial  pumping’  (Price,  1983).  The  area  shown  in  Figure  7 
covers  a  little  less  than  one  day  of  time  along  the  track,  and  so  we  see  only  the  first 
cycle  of  the  inertial  pumping.  Because  Norbert  was  intense  and  comparatively  small, 
the  inertial  pumping  had  a  large  amplitude,  the  simulations  suggesting  that  r\  w  40  m 
at  maximum,  or  more  than  half  the  SML  thickness.  The  largest  upwelling  seen  in  the 
data  was  about  25  m,  but  the  simulations  suggest  that  larger  upwelling  (and  currents) 
may  have  occurred  just  behind  the  hurricane  in  an  area  not  sampled  by  AXCPs. 

The  dynamical  effect  of  the  inertial  pumping  is  to  produce  a  hydrostatic  pressure 
perturbation  within  the  SML  and  the  main  thermocline.  The  region  of  upwelling  just 
behind  the  hurricane  center  causes  a  low  pressure  perturbation  whose  amplitude  in 


Figure  7.  Upwelling  (top  two  panels)  and  thermocline-depth  current  (bottom  two  panels) 
from  Norbert.  The  maximum  upwelling  in  the  simulation  (top  right  panel)  occurs  just 
between  the  AXCP  sections  made  behind  the  hurricane. 
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the  SML  or  upper  thermocline  is  approximately  g'  77  p,  and  decreasing  with  depth  in 
the  thermocline  (Table  III  shows  how  77  depends  upon  external  parameters).  The 
horizontal  scale  of  77  is  roughly  3  Rmax  (Figure  7),  and  thus  the  thermocline  is  accelerated 
toward  the  upwelling  maximum  which  occurs  just  behind  the  hurricane  center.  The 
amplitude  of  the  acceleration  is  A  ss  «  2.5  x  10-5  m  s-2  (see  Fig.  7  of  Price 

(1983)  for  the  horizontal  structure).  This  acceleration  oscillates  in  time  along  with  the 
inertial  pumping,  and  so  the  resulting  thermocline  currents  also  have  a  near-inertial 
time  dependence,  which  is  apparent  in  Figure  7  (lower  panels)  as  a  clockwise  turning 
with  increasing  distance  behind  the  hurricane  center.  The  maximum  amplitude  of 
the  thermocline  currents  is  roughly  the  time  integrated  acceleration  over  the  first  half 
inertial  period,  or,  A  x  2  x  104  s  ss  0.5  ms-1,  which  is  an  appreciable  fraction,  «  0.4, 
of  the  SML  current  (this  fraction  being  proportional  to  and  nearly  equal  to  the  SML 
Burger  number).  Note  that  the  thermocline-depth  currents  have  a  rather  complex  and 
small  scale  horizontal  structure  compared  with  the  very  simple  pattern  of  either  the 
SML  current  or  of  the  upwelling  itself,  since  they  are  forced  by  the  horizontal  gradient 
of  the  pressure. 

The  vertical  structure  of  the  thermocline  response  is  also  very  simple  during  at 
least  the  early  part  of  the  relaxation  stage  that  we  see  here.  The  inertial  pumping 
has  a  nearly  uniform  phase  and  amplitude  through  the  thermocline;  the  pressure  per¬ 
turbation  has  the  same  phase  but  decreasing  amplitude  with  increasing  depth.  The 
thermocline  currents  forced  by  this  pressure  perturbation  should  have  just  the  same 
vertical  structure,  and  clear  evidence  of  this  appears  in  Norbert  AXCPs  13  and  20, 
which  have  the  largest  thermocline  currents  (Figure  8).  They  show  a  very  simple 
structure  in  the  thermocline,  with  nearly  uniform  phase  and  decreasing  amplitude  with 
depth.  Note  though  that  the  phase  difference  between  the  SML  current  and  the  upper 
thermocline  current  can  be  180  degrees,  as  in  AXCP  N13,  or  nearly  zero  as  in  N21. 
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AXCP  N21 


AXCP  N13 


Figure  8.  Current  and  temperature  profiles  from  Norbert  AXCP  21  (upper  left)  and 
AXCP  13  (upper  right).  These  profiles  were  taken  to  the  left  and  right  of  the  track  (see 
Figure  3).  In  this  figure  only,  the  AXCP  data  are  shown  in  a  nearly  raw  form;  the  only 
data  processing  being  a  simple  boxcar  smoothing  over  a  roughly  10  m  interval.  The 
surface  wave  component  has  not  been  removed,  but  was  fairly  small  in  both  profiles.  The 
bottom  two  panels  are  the  model-simulated  current  and  temperature  at  the  locations  of 
AXCP  20  and  13.  Note  that  the  upper  thermocline  current  (depths  below  70  m)  has  an 
appreciable  amplitude,  and  only  small  phase  change  with  depth. 
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This  very  simple  vertical  structure  can  be  simulated  by  a  model  having  only  a  few 
active  thermocline  layers  (compare  the  simulated  current  profiles  in  Figure  8  (lower) 
with  the  observed  profiles  above).  Indeed,  a  model  having  only  one  active  thermo¬ 
cline  layer  can  simulate  the  vertical  phase,  near-inertial  time-dependence,  and  the 
thermocline-averaged  amplitude  of  the  initial  relaxation  stage  response.  However,  as 
the  relaxation  stage  response  continues  to  develop  during  the  next  several  days  to 
weeks,  this  structure  will  become  somewhat  more  intricate  (Price,  1983).  Brink  (1989) 
has  observed  the  thermocline  depth  currents  in  the  Gloria  case  in  moored  array  data, 
and  noted  that  the  phase  change  across  the  thermocline  grew  to  more  than  half  a  cycle 
within  about  a  week  after  the  hurricane  passage.  To  simulate  this  structure  without 
serious  error  requires  more  vertical  resolution  than  is  given  by  one  or  a  few  layers. 

In  the  Norbert  case  the  observed  thermocline  currents  were  up  to  0.3  m  s_1  di¬ 
rectly  under  the  hurricane,  and  model  results  suggest  that  stronger  currents  may  have 
occurred  in  a  region  behind  the  hurricane  that  was  not  sampled  by  AXCPs.  The  en¬ 
ergy  and  momentum  of  these  thermocline  currents  comes  from  the  wind-driven  SML 
current  by  the  mechanism  of  inertial  pumping  and  pressure-coupling  described  above, 
and  this  process  causes  the  SML  current  to  decay  in  time  with  an  e-folding  scale  of 
typically  5  days  (Price  (1983);  this  is  the  process  often  parameterized  in  1-D  models 
by  a  linear  drag,  see  for  example  Pollard  and  Millard  (1970)  who  first  described  this). 
If  the  pressure  coupling  is  arbitrarily  suppressed  in  the  3-D  model,  then  the  simulated 
maximum  SML  current  in  the  Norbert  case  changes  considerably,  increasing  by  about 
35%.  Thus  the  non-local  dynamics  of  the  relaxation  stage  response  were  important 
even  during  the  forced  stage  response  to  Norbert  (or  more  to  the  point,  in  a  case  with 
Burger  number,  M  >  |). 
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7  Vertical  Shear  Within  the  Surface  Mixed-Layer 


In  this  section  we  examine  the  vertical  shear  of  the  SML  current,  and  compare 
AXCP  profile  measurements  with  simulations  of  the  current  made  by  the  1-D  model.  In 
Part  I  we  discussed  at  length  the  problem  of  separating  the  observed  current  measured 
by  AXCPs  into  surface  wave,  shear  and  layer  average  mean  currents.  In  brief,  we  found 
that  the  fitting  procedure  used  to  analyze  the  profiles  gives  very  reliable  estimates  of 
the  layer- averaged  mean  current;  sensitivity  studies  showed  that  bias  of  the  estimated 
mean  SML  current  was  0.02  m  s_1  and  that  the  coefficient  of  variation  was  0.07.  The 
vertical  shear  of  the  SML  current  was  found  to  be  more  difficult  to  estimate  because 
the  change  in  wave  amplitude  with  depth  can  be  aliased  to  appear  as  shear,  especially 
when  the  SML  is  shallow  and  only  one  or  fewer  wave  cycles  are  sampled  within  the 
SML.  The  coefficient  of  variation  for  the  estimated  shear  in  the  SML  was  thus  very 
high,  1.02,  but  the  bias  was  only  —0.0004.  In  comparing  the  shear  estimated  from  the 
AXCP  measurements  with  the  simulations  from  the  1-D  model,  we  can  thus  expect  that 
the  scatter  of  individual  estimates  will  be  very  large  even  if  the  model  were  perfect,  but 
that  the  mean  value  and  mean  trend  over  an  ensemble  of  estimates  should  be  useful 
for  verifying  the  model. 

Figure  9  shows  the  measured  and  simulated  profiles  of  current  magnitude  for 
AXCP  27  from  Gloria.  The  layer  averaged  SML  current  in  this  profile  was  1.70  m  s-1, 
which  was  the  largest  SML  current  measured  in  the  three  storms.  The  shear  in  the 
SML  was  more  or  less  aligned  with  the  mean  current,  and  was  about  0.01  s~J.  Hence, 
the  estimated  surface  current  in  this  profile  was  2.04  m  s-1.  The  shear  in  this  profile  is 
thus  a  significant  feature  for  many  design  studies.  Note  that  the  shear  in  the  transition 
layer  is  only  slightly  larger  than  the  shear  in  the  mixed  layer  for  this  profile.  More 
typically,  the  shear  in  the  stratified  transition  layer  is  considerably  larger  than  is  the 
shear  in  the  SML. 
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Figure  9.  Profile  of  current  speed  measured  at  Norbert  AXCP  27  (solid  line),  and  as 
simulated  by  the  one-dimensional  boundary  layer  model  (dashed  line). 
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The  current  profile  from  the  1-D  model  simulation  is  shown  as  the  dashed  line  in 
Figure  9.  Within  the  SML  and  transition  layers  the  simulated  shear  is  very  similar  to 
that  observed  (in  the  1-D  model  there  is  no  current  below  the  transition  layer  since 
there  is  no  pressure  coupling  with  the  thermocline.)  In  this  and  in  many  cases,  the 
overall  depth  of  vertical  mixing,  as  judged  by  the  thickness  of  the  wind-driven  layer,  is 
reasonably  consistent  between  the  1-D  model  simulation  and  the  observations. 

A  scatter  plot  of  simulated  vs.  measured  SML  shear  in  the  direction  of  the  local, 
flight-level  wind  is  given  by  Figure  10,  which  includes  all  three  cases.  The  simulated 
shear  was  estimated  from  the  1-D  model  results  by  a  fitting  procedure  much  like  that  ap¬ 
plied  to  the  AXCP  data.  The  solid  circles  indicate  flight  level  winds  less  than  30  m  s-1, 
and  the  open  circles  indicate  weaker  winds.  The  scatter  in  the  plotted  points  is  very 
large,  but  expected  from  the  uncertainty  in  the  fit  to  the  measured  profiles  (the  coeffi¬ 
cient  of  variation  between  the  measured  and  simulated  shears  is  about  1,  or  comparable 
to  the  coefficient  of  variation  found  in  the  sensitivity  studies  of  Part  I).  The  same  data 
are  also  shown  in  Figure  11,  where  measured  shear  values  are  open  circles,  and  simu¬ 
lated  shear  values  are  asterisks.  The  solid  and  dashed  lines  are  averages  over  10  m  s_1 
ranges  of  wind  speed  for  the  measured  and  simulated  shears  respectively.  Note  that 
there  are  some  estimates  from  both  the  data  and  the  model  indicating  a  shear  that  is 
opposed  to  the  local  wind.  These  arise  from  cases  where  there  has  been  a  rapid  change 
in  wind  direction,  and  indicate  that  the  shear  and  the  wind  stress  are  not  in  a  steady 
state. 

Despite  the  scatter  and  the  nonstationary  effects,  there  is  a  clear  trend  for  shear 
to  increase  with  wind  speed  in  both  the  observations  and  the  1-D  model  results.  At 
moderate  wind  speeds,  <  20  m  s'1,  the  shear  is  roughly  0.03-1,  and  is  somewhat  over¬ 
estimated  by  the  model.  At  the  largest  wind  speeds,  up  to  40  m  s_1  which  corresponds 
to  a  stress  of  about  4  Pa,  the  measured  shear  is  roughly  0.01  s-1,  and  is  simulated  well 
by  the  1-D  model.  A  shear  of  this  magnitude  can  be  significant  for  some  purposes, 
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Figure  10.  Vertical  shear  of  the  SML  current  in  the  direction  of  the  local  wind.  1-D 
simulated  values  are  plotted  against  measured  values.  The  open  circles  indicate  estimates 
made  where  wind  speed  was  greater  than  than  25  m  s_l,  solid  circles  are  estimates  at  all 
lower  wind  speeds. 
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Figure  11.  Vertical  shear  of  the  SML  current  as  a  function  of  wind  speed.  Open  circles  are 
estimates  made  from  AXCP  data,  and  asterisks  are  estimates  made  from  the  1-D  model 
simulations.  Solid  and  dashed  lines  are  average  trends  computed  from  the  observations 
and  simulations  respectively. 


and  specifically,  it  shows  that  a  correct  layered  model  would  tend  to  underestimate  the 
surface  currents  by  about  0.2  ms-1  in  hurricane  conditions. 


8  Summary  and  Remarks 

We  believe  that  the  most  important  result  of  this  study  and  Part  I  is  the  demon¬ 
stration  that  useful  observations  of  ocean  currents  can  now  be  acquired  under  the  most 
extreme  weather  conditions,  and  that  numerical  ocean  models  are  able  to  rationalize 
the  main  features  of  such  data  sets.  Our  practical  goal,  to  observe  and  model  storm- 
driven  currents  for  design  purposes,  appears  to  be  within  reach,  at  least  for  open  ocean 
conditions. 

The  scientific  goals  for  this  work  were  somewhat  more  open-ended,  but  we  hope 
partially  realized  here  by  showing  how  the  horizontal  structure  of  the  forced  stage  SML 
response  is  imposed  by  the  atmospheric  forcing  due  to  the  hurricanes.  If  a  storm  has  a 
non-dimensional  speed  S  of  0(1),  as  these  and  most  hurricanes  do,  then  the  dominant 
horizontal  structure  of  SML  currents  can  be  understood  as  the  local  response  of  the 
ocean  to  a  time  variable  wind  stress.  At  locations  where  the  stress  turns  in  the  sense  of 
an  inertial  current  (right  side  of  the  track),  the  response  can  be  greatly  enhanced.  The 
across-track  scale  of  the  response  is  the  scale  of  the  hurricane,  and  the  along-track  scale 
is  just  the  product  of  the  storm  translation  speed  and  the  inertial  period.  The  vertical 
structure  of  the  forced  stage  response  is  also  quite  simple.  The  wind-driven  SML  current 
has  a  fairly  small  vertical  shear,  typically  0.01  s-1  under  very  high  stress  conditions, 
which  may  nonetheless  be  important  for  some  design  purposes.  Most  vertical  shear 
occurs  within  a  stratified  transition  layer  which  has  a  thickness  roughly  half  that  of  the 
SML. 
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In  the  Norbert  case,  which  is  characterized  by  a  fairly  large  SML  Burger  num¬ 
ber,  we  can  see  the  start  of  the  relaxation  stage  response;  appreciable  upwelling  and 
thermocline-depth  currents  occurred  beneath  the  trailing  edge  of  the  hurricane.  These 
thermocline-depth  currents  also  have  a  simple  vertical  structure,  having  very  little 
phase  change  with  depth  and  a  monotonically  decreasing  amplitude.  Their  horizontal 
structure  is  more  complex,  but  is  directly  related  to  the  horizontal  structure  of  the 
forced  stage  response  through  the  field  of  inertial  pumping  and  hydrostatic  pressure 
perturbation. 

Several  specific  modeling  questions  were  raised  in  this  study,  and  were  addressed 
by  comparison  of  model  simulations  with  the  observations.  First,  we  found  that  for  the 
purpose  of  modeling  ocean  currents  the  hurricane  wind  stress  fields  can  be  estimated 
well  from  conventional  methods.  Sub-hurricane  scale  errors  may  be  present  (in  the  left 
rear  quadrant),  but  are  not  readily  identified  in  these  data  sets.  Second,  we  found  that 
the  bulk  Richardson  number  estimated  from  AXCP  data  takes  on  values  just  less  than  1 
in  regions  near  the  hurricane  center  where  vertical  mixing  was  expected  to  be  strongest, 
and  is  larger  in  outlying  regions.  This  is  consistent  with  a  Richardson  number  closure 
for  mixing.  Third,  we  noted  here  and  in  Part  I  that  the  forced  stage  response  and 
the  early  relaxation  stage  response  have  a  very  simple  vertical  structure  that  can  be 
represented  by  a  layered  model.  If  the  Burger  number  is  appreciable  (as  in  the  Norbert 
case),  a  3-D  model  having  the  dynamics  of  inertial  pumping  seems  necessary,  even  for 
the  forced  stage  response.  Fourth,  the  high  resolution  1-D  model  was  found  to  give 
useful  simulations  of  the  SML  current  vertical  shear,  especially  in  the  range  of  very 
large  wind  speed  where  the  vertical  shear  is  largest  and  most  important. 

These  data  sets  show  the  great  power  that  aircraft-deployed  expendable  instru¬ 
ments  can  bring  to  studies  of  air/sea  interaction.  Aircraft  can  provide  oceanographic 
sampling  in  definite  relationship  to  synoptic  weather,  and  at  the  same  time  provide 
a  high  quality  view  of  the  winds.  However,  our  problems  with  interpretation  of  the 
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Sargasso  Sea  data  sets  should  be  a  caution  to  future  investigators  that  one-time  sur¬ 
vey  data  will  not  always  be  sufficient  to  sort  out  storm-induced  currents  from  ambient 
currents.  In  future  studies  it  would  be  useful  to  acquire  initial  survey  data,  and  just  as 
important  and  perhaps  more  effective,  avoid  regions  of  high  ambient  variability.  The 
AXCP  sampling  plan  employed  in  these  studies  appears  to  have  sampled  the  SML 
currents  fairly  well,  especially  the  regions  of  strongest  current.  Sampling  within  the 
thermocline  appears  to  have  been  less  satisfactory,  and  could  be  easily  improved  in 
future  studies  by  making  a  section  across  the  hurricane  track  and  about  100  km  behind 
the  hurricane  center. 
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Appendix 

AXCP  Data  Tabulations  and  Analysis  Method 

a)  AXCP  Station  Data 

The  P3  aircraft  flew  a  similar,  more  or  less  star-shaped  pattern  through 
each  hurricane  (Fig.  1A  shows  the  Gloria  flight  path).  AXCPs  were  dropped  at 
roughly  25  km  intervals,  with  enhanced  sampling  on  the  right  side  of  the  track.  In 
each  case  about  15  AXCPs  produced  useful  data  (Table  1A). 

The  AXCPs  measure  temperature,  and  the  motionally-induced  electric  field 
set  up  by  ocean  currents  and  the  orbital  motions  of  surface  gravity  waves.  They 
fall  through  the  water  at  a  rate  W  =  4.5  m  s-1  to  a  depth  of  1500  m  (fast-fall 
probes),  or  at  a  reduced  rate,  W  =  2.2  m  s_1,  to  a  depth  of  200  m  (slow-fall 
probes).  Because  of  a  rapid  wind-driven  drift  of  the  surface  buoy  due  to  very  high 
winds  encountered  in  hurricanes,  most  of  the  fast-fall  probes  produced  usable  data 
down  to  only  about  800  m,  which  is  more  than  adequate  for  this  study. 

b)  AXCP  Analysis 

The  currents  inferred  directly  from  AXCP  measurement  are  relative 
currents  because  of  an  unknown  reference,  f/r,  which  is  independent  of  depth.  In 
these  experiments  we  have  estimated  Ur  as  the  depth-independent  current  in  the 
deepest  portion  of  the  profiles  (see  Figure  2A  for  an  example).  The  estimated 
reference  currents  were  usually  less  than  10  cm  s-1,  and  their  uncertainty  is  not 
thought  to  be  significant. 
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Table  1A:  AXCP  station  data. 


D  #  is  the  AXCP  number;  N  denotes  Norbert,  J,  Josephine.  G,  Gloria.  F/S  indicates  fast-  or  slow-fall  probe  V  and  Y 
are  the  along-track  and  across-track  coordinates  in  the  storm-centered  system;  X  increases  toward  the  top  of  the  figure 
V  increases  to  the  right  TML  is  temperature  of  the  surface  mixed  layer,  and  ST  is  the  temperature  change  across  the 
transition  layer,  r,  is  the  apparent  upweiling  within  the  thermocline.  Z,  is  the  depth  of  layer  i 
a:  AXCP  G19  did  not  go  deep  enough  to  measure  layer  3. 


D  #  F/S  lat. 

deg.  N 


N  2 

F 

20  56 

N  3 

F 

18.11 

N  4 

F 

20.25 

N  6 

F 

19.60 

N 13 

F 

19.07 

N 14 

F 

18.96 

N 15 

S 

19.41 

N 16 

S 

19.20 

N 18 

F 

18.83 

N20 

F 

18.51 

N21 

F 

18.33 

N22 

S 

18.88 

N23 

s 

18.90 

N24 

s 

18.91 

N26 

F 

18.93 

N3I 

F 

19.81 

J  2 

F 

28.61 

J  3 

F 

28.83 

J  4 

F 

29.03 

J  7 

F 

29.37 

J  8 

S 

29.07 

J13 

F 

29.28 

J 14 

S 

29.30 

J17 

s 

29.33 

J20 

s 

30.40 

J21 

F 

30.66 

J25 

F 

27.93 

J26 

F 

27.75 

J27 

F 

27.57 

J29 

F 

27.41 

G  5 

S 

28.31 

G  7 

s 

28.15 

Gil 

s 

29.86 

G12 

s 

29.40 

G13 

s 

28.46 

G15 

s 

28.13 

G16 

s 

30.23 

G17 

F 

30.43 

G18 

F 

29.95 

G19 

F 

29.79 

G21  S  29.47 

G24  S  28.93 

G25  F  28.91 

C27  S  28.87 

G31  F  28.68 


Ion.  date 

deg.  E 

-108.20  23SEP84 

-108.76  24SEP84 

-107.86  23SEP84 

-108.33  23SEP84 

-108.31  23SEP84 

-108.07  23SEP84 

-109  08  24SEP84 

-108.95  23SEP84 


-108.88 

24SEP84 

-108.83 

24SEP84 

-108.80 

24SEP84 

-108.93 

24SEP84 

-108.83 

24SEP84 

-108.93 

24SEP84 

-109.68 

24SEP84 

-109.40 

24SEP84 

-73.77 

110CT84 

-73.48 

110CT84 

-73.20 

110CT84 

-72.28 

110CT84 

-72.22 

UOCT84 

-70.52 

110CT84 

-70.70 

110CT84 

-71.97 

110CT84 

-71.38 

110CT84 

-71.20 

110CT84 

-71.76 

UOCT84 

-71.65 

110CT84 

-71.55 

UOCT84 

-72.23 

110CT84 

-74.09 

26SEP85 

-73.90 

26SEP85 

-74.22 

26SEP85 

-74.78 

26SEP85 

-75.87 

26SEP85 

-76.26 

26SEP85 

-73.83 

26SEP85 

-73.61 

26SEP85 

-74.89 

26SEP85 

-74.89 

26SEP85 

-74.91 

26SEP85 

-74.92 

26SEP85 

-73.16 

26SEP85 

-73.85 

26SEP85 

-77.00 

26SEP85 

time 

X 

Y 

UTC' 

km 

km 

2244 

149.2 

65.8 

17 

-61.5 

-113.7 

2253 

157.5 

14.4 

2311 

74.0 

-8.7 

2335 

40.6 

-58.5 

2339 

53.7 

-85.5 

112 

0.0 

0.0 

2359 

-4.6 

-8.3 

5 

-23.5 

-44.4 

11 

-40.7 

-75.2 

13 

-50.2 

-92.2 

30 

-23.8 

-43.0 

31 

-14.0 

-48.8 

31 

-21.8 

-40.8 

101 

-83.6 

7.0 

123 

-0.1 

52.3 

1306 

-130.5 

-147.7 

1300 

-103.1 

-119.8 

1255 

-76.3 

-94.2 

936 

0.0 

0.0 

931 

111 

-26.3 

1004 

196.7 

19.2 

1001 

176.4 

18.3 

943 

35.3 

0.6 

1041 

87.2 

101.2 

1054 

103.8 

126.6 

912 

78.9 

-122.0 

909 

93.7 

-136.4 

903 

107.1 

-150.3 

854 

34.5 

-176.7 

950 

13.7 

-125.0 

946 

26.8 

-147.6 

833 

95.2 

34.6 

842 

16.4 

20.8 

903 

-140.1 

-10.9 

910 

-195.7 

-21.6 

825 

152.8 

49.0 

821 

184.8 

56.2 

604 

82.7 

112.1 

620 

70.3 

94.4 

558 

61.5 

72.9 

550 

39.3 

28.4 

729 

179.5 

-85.9 

720 

112.5 

-52.4 

637 

-193.4 

100.4 

TML 

ST 

n 

C 

C 

m 

28.4 

2.5 

-5. 

27.2 

5.0 

20. 

28.5 

3.0 

0. 

27.8 

3.0 

0. 
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0. 

26.5 

3.0 

15. 

27.4 
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27.5 

50 

-5. 

27.0 

5.0 

5. 

26.7 

5.0 

5. 

27.2 

8.0 
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4.5 
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Figure  2A.  Three-layer  model  fit  to  an  AXCP  current  and  temperature  profile. 
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For  most  of  our  purposes  it  was  necessary  to  separate  the  currents  from  the 
wave  orbital  motions,  and  it  was  convenient  to  work  with  layer- averaged  current 
and  shear  rather  than  an  arbitrary  profile.  We  have  therefore  fitted  the  observed 
profiles  to  a  model  comprised  of  a  single  surface  wave  and  a  three-layer  linear 
current  profile,  L(z).  For  the  east  component, 

U(z)  =  ekz  [a  cos(uz/W)  +  6sin(u>z/W)]  +  L(z)  , 

where  k  =  uj2 / g  is  the  vertical  wavenumber,  u>  is  the  wave  frequency  estimated 
from  the  observed  wave  oscillation  and  given  the  known  fall  rate,  W ,  and  g  is  the 
acceleration  due  to  gravity.  The  coefficients  a,  b  are  found  by  the  fitting,  and  note 
that  the  surface  amplitude  of  the  wave  is  just  JJw  =  (a2  +  b 2)1^2.  L(z )  is  the 
current  profile 

=  u,  +  Uzi  (^  -  f-)  if  Z0  >  z  >  Z1 

L(z )  =  U2  +  UZ2  (Z  -  Zl+Z 2)  ^  Zi  >  z  >  Z2 

=  u3  +  UZ3  [z  -  if  Z2>z>Z3 

where  Zq  is  the  start  of  usable  data,  Z\,  Z2,  Z3  are  the  depths  of  the  layers,  and 
U,  and  Uzi  are  the  depth-averaged  current  and  shear  found  in  layer  i  by  the  fitting 
procedure.  They  are  constrained  to  yield  a  continuous  profile,  i.e.  at  the  base  of 
layer  1, 

Ul  +  £i^i  =  Ul  +  u„^zIi  . 

The  layer  depths  Z,  were  chosen  subjectively  based  upon  the  observed 
structure  of  the  temperature  and  current  profiles,  and  with  a  definite  physical 
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model  in  mind.  Layer  1  is  the  surface  mixed  layer  (SML)  over  which  temperature 
is  uniform  though  current  may  not  be;  layer  2  is  the  transition  layer  at  the  base  of 
the  SML  which  is  strongly  sheared  and  stratified;  layer  3  is  the  upper  thermocline. 
This  three-layer  model  seemed  apt  for  most  of  the  profiles,  and  especially  those 
with  strong  hurricane-driven  currents.  This  can  be  judged  in  part  by  noting  that 
the  root  mean  square  current  which  could  not  be  accounted  by  the  fit,  R,  was 
typically  about  8  cm  s_1  (see  Tables  1  and  2  of  Part  1,  and  Table  2A  here). 

The  only  troublesome  aspect  of  the  fitting  procedure  arose  in  cases  where 
we  had  a  fast-fall  AXCP  and  a  shallow  SML  (most  common  in  Norbert).  In  that 
circumstance  the  shear  within  the  SML  and  the  wave  velocity  may  be  almost 
indistinguishable.  Because  they  are  not  orthogonal,  the  fitting  procedure  has  a 
tendency  to  return  compensating  estimates  for  shear  and  wave  coefficients  which 
are  likely  to  be  too  large.  To  investigate  this,  we  carried  out  a  series  of  sensitivity 
tests  using  synthetic  current  profiles  (reported  in  Part  1).  We  found  that  the 
coefficient  of  variation  for  shear  estimates  is  large,  about  1,  but  that  the  estimates 
were  unbiased.  We  expect  a  simlar  result  to  obtain  in  the  present  analysis; 
individual  shear  estimates  may  have  0(1)  uncertainty,  but  we  expect  that  an 
average  over  many  samples  will  be  unbiased. 

c)  Storm- Centered,  Coordinate  System 

The  AXCPs  were  dropped  over  a  period  of  about  four  hours  during  which 
the  hurricanes  moved  a  significant  distance.  In  order  to  produce  a  quasi-synoptic 
field  we  have  therefore  plotted  AXCP  data  in  a  storm-centered  coordinate  system 
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Table  2A:  Coefficients  for  three-layer  model  fit  to  Gloria  AXCP  profiles. 

D  #  is  the  AXCP  number.  Z0  is  the  start  of  usable  data  (m)  TML  is  the  temperature  of  the  surface  mixed  layer  (C) 
P  is  the  period  of  the  surface  wave  (s),  and  a  and  b  are  the  cosine  and  sine  coefficients  of  the  fit  to  the 
surface  wave  (cm  s'1).  The  first  row  of  data  for  each  AXCP  is  the  east  component  and  the  second  row  is  the 
north  component  (all  directions  are  true).  Z,  is  the  depth  of  the  bottom  of  layer  i  (m), 

I  \  is  the  mean  current  in  the  layer  (cm  s'1),  and  U,x  is  the  shear  (cm  s'1  m'1).  The  reference  velocity  has  been 
subtracted  away  R  is  the  rms  difference  between  the  observed  and  the  best  fit  profile  (cm  s'1). 

The  format  of  this  table  exactly  follows  Tables  1  and  2  of  Part  1  which  give  Norbert  and  Josephine  data. 
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whose  origin  is  the  hurricane  eye  position  Xc  at  the  center  time  of  the  survey,  tc 
(Table  II).  Thus  a  drop  made  at  time  tj  =  tc  +  A t  and  at  position  Xd  —  Xc  +  A.Y 
would  have  a  storm-centered  coordinate 

X  =  Xd  -  Xc  +  A*  Uh  , 

where  Uh  is  the  hurricane  translation  speed  (assumed  constant  during  the 
duration  of  the  flight). 

Finally,  to  simplify  comparison  of  the  three  cases,  the  coordinate  X  is 
rotated  into  a  frame  in  which  Uh  is  up  the  page.  The  coordinates  are  then  termed 
along-track  and  across-track,  and  listed  in  Table  1A  as  (X,  Y)\  the  corresponding 
current  and  shear  components  are  in  Table  3A. 
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Table  3A:  AXCP  current  data  in  tbe  storm-centered  coordinate  system. 

D  #  is  the  AXCP  number;  N  denotes  Norbert,  J,  Josephine,  G,  Gloria.  Uw  and  Vw  are  the  along-  and  across-track 
amplitude  of  the  wave  component  at  the  surface  (direction  is  ambiguous  to  180°).  Ui  and  V,  are  the  along-  and  across-track 
currents  averaged  over  layer  i,  and  U„,  are  the  corresponding  shears, 
a:  AXCP  G19  did  not  go  deep  enough  to  measure  layer  3. 
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